function [x,k]=SORmethod(A,b,x0,N,emg,w)
%A是线性方程组的左端矩阵
%b是右端向量
%x0是迭代初始值
%N表迭代次数上限，若迭代次数大于N，则迭代失败% emg 表示控制精度
% w表示松弛因子
%用SOR迭代法求线性方程组Ax=b的解
%k表迭代次数
%x表示用迭代法求得的线性方程组的近似解n=length(A);
x1=zeros(n,1)；x2=zeros(n，1)；
xl=x0；
r=max(abs(b-Axl))；
k=0； 附录0MATL
while r＞emg
for i=1:n
sum=0；
for j==1:n
if j>=i
um=sum＋A(i,j)*x1(j)；
elseif j<i
sum=sum+A(i,j)*x2(j)；
end
end
x2(i)=x](i)+w×(b(i)-sum)/A(i，i)
end
r-max(abs(×2-x1))；
u 
x1=x2;
k-k+1;
if k>N
disp(迭代失败，返回)；
rcturn;
end
end
x=x1；